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Abstract 



We introduce a theory for the absorption of electromagnetic radiation by 
small metal particles, which generalises the random phase approximation by 
incorporating both electric and magnetic dipole absorption within a unified 
self-consistent scheme. We demonstrate the equivalence of the new approach 
to a superficially dissimilar perturbative approach. We show how to obtain 
solutions to the self-consistent equations using a classical approximation, tak- 
ing into account the non-locality of the polarisability and the conductivity 
tensor. We discuss the nature of the self-consistent solutions for diffusive and 
ballistic electron dynamics. 
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I. INTRODUCTION 



The quantum theory of absorption of radiation by small metal particles has proved to 
be a difficult area, partly because it has not always been realised that the internal electric 
field must be treated self-consistently, and partly because of confusion about how to obtain 
the self-consistent field. The literature is surveyed in [1-3], and the introductory sections of 
[4] and [5] include a brief review of the literature relevant to the present paper. This paper 
extends a series of works [4-8] which have considered the self-consistent potential from a 
semiclassical viewpoint. 

First we clarify the regimes in which our results are applicable. The following frequency 
scales are relevant to the interaction of small metal particles with radiation. The lowest scale 
is given by the mean level spacing A, i.e., wa = A/h. In d dimensions oc a~ d , where a is 
the characteristic size of the particle. Another typical frequency scale is given by the inverse 
of the typical time taken for an electron to traverse the particle. In systems with diffusive 
electron motion it is given by u c = D/a 2 , where D is the diffusion constant. In particles with 
ballistic electron motion, oj c is given by the inverse time of flight, oj c = vp/a, where vp is the 
Fermi velocity. The highest frequency scales are the plasma frequency uj p and the frequency 
ljf = Ep/h derived from the Fermi energy. In two dimensions, u p oc aT 1 / 2 , whereas in three 
dimensions u v is comparable to the Fermi frequency, which is independent of the size of the 
particle. The frequency scales are therefore ordered as follows 

<^ <^ < • (1) 

Our results are relevant to frequencies satisfying u ^> c^a (which justifies the neglect of 
quantum effects), and w<w p (which means that screening of external fields by polarisation 
of the particle is significant). There are six relevant length scales in the problem, namely 
the wave length A of the external radiation, the linear dimension a of the particle, the skin 
depth A s , the Thomas- Fermi screening length A T f, the Fermi wavelength A F and the mean 
free path /. In the following it will be assumed that 

A, A s > a > A T f ~ A F . (2) 

Both diffusive (I <C a) and ballistic (I ^> a) dynamics will be discussed. 

Under these conditions the absorption can be divided into electric and magnetic dipole 
contributions. In earlier papers these were discussed separately: the electric dipole absorp- 
tion was discussed for various different situations in [4-7], and the magnetic dipole absorption 
case was discussed in [8]. It is desirable to have a unified description. In section II we de- 
velop a generalisation of the commonly used 'random phase approximation' [9] scheme which 
gives a uniform approach to both electric and magnetic dipole absorption. Most treatments 
of the electric dipole absorption coefficient are based upon calculating the imaginary part 
of the polarisability. References [4-7] used an alternative approach, namely time dependent 
perturbation theory in the effective potential. In section III we show that these apparently 
dissimilar approaches are equivalent. 

We also discuss a unified approach to determining the self-consistent field: we present 
some new results for the ballistic case, where the bulk mean free path of the electrons is 
large compared to the dimensions of the particle. Section VI describes the form of the non- 
local polarisability, and section V discusses a general semiclassical method for determining 



2 



the self-consistent field. We show that the self-consistent potential used in [4,6], which 
treated electric dipole absorption ballistic systems, was not correct for u> ^> oo c . Section VI 
summarises the results, and presents an argument indicating that the results of [4,6] are 
nevertheless qualitatively correct. 

II. CALCULATION OF THE ABSORPTION COEFFICIENT 

This section will discuss the general principles underlying the calculation of the electro- 
magnetic response. 

A. Formulation of the problem 

An electromagnetic wave induces currents which result in both electric and magnetic 
polarisation of a conducting conducting particle. In what follows we will only consider 
linear effects (where the polarisation is proportional to the applied field), and the externally 
applied field will be assumed to be uniform over the dimension of the particle. We will only 
be concerned with the coefficients relating dipolar moments to the externally applied field: 
higher moments will not be considered. The electric dipole d and magnetic dipole m of a 
single particle are given by 

d = a E cxt , m = f3 B cxt + 7 E cxt (3) 

where E ext and B cxt are the externally applied electric and magnetic fields, a and (5 are the 
electric and magnetic susceptibility tensors of the particle. The cross-susceptibility 7 is not 
usually included. It is absent for spherical particles and some other symmetric geometries, 
and when it is non-zero it vanishes in the low- frequency limit. We will ultimately give a 
completely general treatment, showing that the cross-susceptibility makes no contribution 
to the absorption, but for simplicity the cross-term will be dropped in the remainder of this 
introductory section. The externally applied fields are assumed to be multiplied by a factor 
of the form exp(— icot), and the polarisability tensors are understood to be functions of 00 
with complex valued components, because there may be a phase shift between the applied 
field and the response. For example, the actual value of the dipole moment at time t is taken 
to be d(t) = Re[dexp(— iut)]. 

These polarisations are detectable at a macroscopic level in various ways: they alter 
the dielectric constant and magnetic permeability of the medium in which the particles 
are dispersed, and they may also be detected by observing scattering and absorption of 
radiation. The polarisability determines two processes which result in the attenuation of 
radiation, namely scattering and absorption of energy. Both of these processes can be 
characterised conveniently at the microscopic level by calculating the rate of loss of energy 
from the incident beam due to interaction with a single particle: this will be denoted by 
dE/dt [where X denotes the time average of X(t)]. The two most commonly used measures 
used to quantify these processes are the cross section per particle S and the attenuation 
coefficient 7. To relate the energy loss to the cross section, note that the energy density in 
an electromagnetic wave is |e E 2 : the cross section is therefore 
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2 dE 

where X denotes the time axerage of X(t). The attenuation coefficient 7 is defined by the 
expression I = J exp(— 72), where I is the intensity at distance z along the beam. The 
attenuation coefficient is given by 7 = AfS, where M is the particle density. 

At low frequencies the real part of the polarisability approaches a constant (and the 
imaginary part approaches zero). It follows that at sufficiently low frequencies the scattering 
cross section scales as u 4 . It will be shown that in the absorption cross section typically scales 
as u 2 , implying that absorption is expected to be the dominant process at low frequencies. 

The absorption of radiation can be related to the imaginary parts of the polarisability 
tensors: we will give a careful explanation of this. Electron spin is not significant in this 
context, and the full Hamiltonian for the electrons is taken to be 

N -1 N N 2 

H = E ^[P* - eA cxt (r,,t)] 2 + V(r t ) + cxt (r,,t) + ± £ £ — — - . (5) 

i=l Zm i=l j=i 47re 0| r i r j\ 

The externally applied electric and magnetic fields are considered to be spatially uniform, 
since the particle is small compared to the wavelength of the radiation, we ignore the spatial 
dependence of the electric and magnetic fields, and write: 

cxt (r,t) = er.E cxt , V A A cxt (r,t) = B cxt (t) . (6) 

In the case where the circularly symmetric gauge 

A cxt = A rot (r, t) = |B cxt (t) A r (7) 

is used, the full Hamiltonian contains terms coupling the system to the electric and magnetic 
fields, of the following form: 

H(t) = H + eX.E cxt (t) + -^L.B cxt (t) + 0(B e 2 xt ) 

2m 

N N 

X = ^r l , L = 5>Ap, (8) 

i=i i=i 

where X and L are the total dipole operator and total angular momentum operators, and 
Hq is independent of time. 

To facilitate the calculations we will consider ensemble averages of quantities: if the 
electron motion is ergodic, this is the microcanonical average, and in general the ensemble 
is defined by the region of phase space explored by the dynamics. Angle brackets will be 
used for the appropriate ensemble average. For a general choice of gauge the instananeous 
rate of absorption is then 

/ dE \ I OH \ I ^ \ f 

\~dt) = \~dt) = e (E v *W-[E C xt(t) + E ind (t)]^ = J dr j(r,t).[E cxt (t) + E ind (r,t)] (9) 
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where Vj is the velocity of the 2 th electron, E in d = i^A cxt is the electric field induced by the 
varying magnetic field, and 

i{T,t) = e(£v i {t)6[T-r i {t)]) (10) 

is the current density within the particle. In the special case where the circularly symmetric 
gauge is used this reduces to: 

(f) = (f) = e<t> - E « (t) + ^ (i> - B » (i) 

where (X) and (L) are suitable averages of the centre of mass and angular momentum 
operators. It is impractical to calculate these averages from the full Hamiltonian (0), and 
in the next sub-section it will be shown how they may be estimated using an effective 
Hamiltonian, containing effective fields A e g and e fr- At this stage we will only assume 
that these averages are proportional to the applied fields. These quantities (X) and (L) are 
related to the electric and magnetic dipole moments d and m: 

d = e(X), m=— (L). (12) 
m 

The rate of absorption is obtained by substituting for the time dependence of a monochro- 
matic field using (Q), and ignoring the cross term: 

(-^-} = — Re[dexp(— hot)] Refiu; E ext exp(— icot)} 

— |Re[mexp(— iut)] Re[iu B ext exp(— iui)] . (13) 

Averaging over time gives the general form for the rate of absorption 



(— ) = ±u,E+ xt (5 - S+)E ext + |u;B+ t (/3 - /3+)B ext . (14) 

In the case where the polarisability tensor is isotropic, and the radiation field is plane 
polarised, this expression becomes 



(~dt) = ^Muuiu)} |Eext| 2 + iu1m\p u (w)] |B cxt | 2 (15) 



We note that under the assumptions listed above, the absorption is expressed as the sum of 
two terms, which are naturally referred to as the electric and magnetic dipole absorption. 
Our final result will not neglect the magnetic dipole moment which may be induced by the 
electric field, but we will show that within the framework of our self-consistent approximation 
scheme the cross term in (|3|) makes no contribution to the absorption. The energy absorbed 
does not accumulate in the system electronic system: most of it is eventually transformed 
into heat by interaction with phonons. 
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B. Self-consistent fields 



The Hamiltonian will be approximated by an effective Hamiltonian, in which the elec- 
trons move independently. The direct interaction with the magnetic field via electron spin 
can also be neglected, and the effective Hamiltonian is of the form 

N 1 

#es = E 7T- [Pi - eA cfr (r,,t)] 2 + Kff(ri) + &ff(*i>*) ■ ( 16 ) 

The response of the system is determined by the interaction of the electrons with the elec- 
tric field inside the particle, which is described by the effective potentials A e g and e fr- The 
effective potentials are themselves determined by the distribution of charge within the par- 
ticle. The external magnetic field is also augmented by an induced magnetic field which is 
produced by the action of the currents which flow in order to establish the electric polarisa- 
tion. Provided the particle is sufficiently small, the induced magnetic field can be neglected, 
and our self-consistent theory will yield and equation for the electric field E(r, to) within the 
particle, which is related to the time-dependent electric field as follows: 

/oo 
d*E(r,t) exp(-icjt) . (17) 
-oo 

This field satisfies the Maxwell equations 

V.E = — , V A E = iwB ext . (18) 

The electric field produced by induction when the external electric field is zero will be 
denoted E inc i. The total effective electric field is 

E(r, u) = E ext (r, to) + E ind (r, u) + V0 po i(r, u) . (19) 

The uniform external electric field satisfies V. E ext = 0, and can be derived from an external 
potential: 

E cxt (r,cu) = V0 ext (r,cu) . (20) 

The potential po i results from polarisation of the particle due to the external electric field, 
and is given by 

P oi (t,u) = / dr — 21 

47reo J |r — r | 

where p po i(r, ui) is the charge density resulting from polarisation induced by the external 
electric field, but excluding any polarisation which may result from the induction field. It 
will be convenient to express (|19|) using the notation 

|E) = V|0 cxt ) + |E ind ) + VU\p pol ) (22) 

where U is an operator defined by (pT|), acting on the 'field vector' |p po i)- The dependence 
upon frequency will usually be shown explicitly for operators, but not for field vectors. 
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The current density j(r, u) flowing in the sample to build up the charge density p(r,u>) 
may be assumed to be linearly related to the electric field E(r, uj) in the sample 

j(r,w) = Jdr'E(r,r';u)E(v',u). (23) 

In condensed notation, we write 

U) = IE) (24) 

where £(u;) is the conductivity operator. The non-local conductivity tensor E(r, r', u) is 
related to the non-local polarisability operator H(r,r',u), which gives the charge density 
induced by a potential <f>(r): we write 

p(r, u) = J dr' n(r, r'; u) <f){r', uj) (25) 

or in condensed notation 

\p) = n(w) |0) . (26) 

The polarisability operator II (o>) can be related to a non-local conductivity tensor £(a;) 
by a continuity relation. Using (|23|) and applying the continuity equation, we find (with 
summation over repeated indices implied) 



= J dr' iu;n(r,r';a/)0(rO - V^Sy (r, r'; o;)V^(r') 



(27) 



assuming that the normal component of £(o>) vanishes on the boundary. Upon integration 
by parts, after noting that the resulting equation is valid for any field <p(r,u), we find 

iu fl(u) = - VS(w) V . (28) 



The quantities II (a>) and S(a;) enable (P2|) to be expressed in terms of the electric field alone, 
yielding a self-consistent equation. We write 

|E) = V|0 cff ) + |E ind ) (29) 

where e ff( r , is an effective potential. We consider the solutions for the field E ind and the 
potential e fr separately. The charge induced by the field E in d is 

bind) = -VS( W ). |E ind ) . (30) 

Applying the first of the Maxwell equations ( |i~8|) and using (|30D gives 

V.[|E fa m)-iwE(w)|E ind )] = (31) 

which is the self-consistent equation which must be solved for the field E ind . For the effective 
potential, we find 

\M = |0ext) + t>nM|0cff) • (32) 

This self-consistent equation is sometimes referred to as the 'random phase approximation' 
[9]. Equations ([31]) and (0) must must be solved for the self-consistent fields. We will 
consider semiclassical methods for solving them in section V. 
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C. The rate of energy absorption 

The rate of energy absorption is given by (|]). Averaging over time gives 



(^) = ! Re f dr j*(r,o;).E ext + iRe jdr j*(r,u;).E ind (r, u) . (33) 
In condensed notation this will be written, by analogy with Dirac notation, as 



Using (H), 



/dE\ 
\~dt/ 



(it) = 2 (j|EeXt) + 2 Re ^ j ' Eind ) • ^ 



±Re(j|E)-±Re(j|Vf>|p pol ) = | Re (j|E) - | Re L;(p pol |f)|p pol ) (35) 



where |p po i) = n|0 e fr); the final equality follows from an integration by parts, and use of 
the continuity equation. Using the fact that U is self-adjoint, we obtain the very simple 
expression for the absorption 



(^)^Re(JIE). (36) 



Using the continuity equation and (|2^) , fl28| ) and (|29|). 



(j|E) = iu (0 cfr |n+H|0 cfr ) + (E ind |S+(^)|E ind ) . (37) 
This gives our final expression for the absorption 



(~dt) = ^ Im (^ff|nM|0 c fr) + |Re(E ind |SH|E ind ) . (38) 

These are two independent contributions to the the rate of absorption, depending on the 
electric and magnetic fields respectively. It is not obvious that these are correctly identified 
as the electric and magnetic dipole coefficients, because the electric field may induce a charge 
density with non-zero angular momentum. We will now show that the first term is due solely 
to the electric dipole. Using the continuity equation and an integration by parts, we find 

Re(j|E ext ) = ^lm(p|0 ext ) = o;Im[d.E ext ] . (39) 

The electrically induced absorption therefore depends only upon the induced dipole moment, 
and is independent of the magnetic moment induced by the electric field. 



III. EQUIVALENCE WITH PERTURBATION THEORY 

A. An alternative expression for the the absorption coefficient 

In this section, we concentrate on the electric absorption. We describe an alternative 
approach to calculating the absorption coefficient, which was used in [4-7], and show that it 
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is equivalent to the first term in fl39|) provided the polarisation operator TL(uj) is related in 
a simple way to a propagator P(uj). This relation will be established in section III B. 

We will consider the action of the effective potential (p e s(r,t) on the electrons. We 
may use either quantum mechanical or classical perturbation theory. We will describe the 
quantum mechanical approach, and will use semiclassical approximations: a classical theory 
in which quantum mechanics only enters in choosing the Fermi-Dirac distribution for the 
initial distribution of electrons gives identical results. Conceptually, the simplest method 
for calculating the absorption in the first using the Fermi golden rule. This is expressed in 
terms of matrix elements cj) nm of the perturbation in the basis \ip n ) formed by the eigenstates 
of the single-particle effective Hamiltonian, H e s'- 

H e $\lp n ) = E n \lp n ) 

<t>nm = (i>n\<i>\i>m), <t> = 4>cs(r) (40) 

The Fermi golden rule states that the rate of transition from an initially occupied state 
to a quasi- continuum of final states, with density of states n and with energy differing by 
Tito from the original state, is given by 

^=^(l0n m | 2 ) (41) 

where the angle brackets denote an average over matrix elements ( < 0n|0eff|V'm)- We will 
consider the case where both the temperature and the photon energy are small compared 
to other energy scales in the problem; generalisations are straightforward. Absorption of 
energy occurs due to the excitation of electrons in occupied states below the Fermi level to 
empty states above the Fermi level. The number of states which can be excited is ~ nhu), 
and the energy absorbed in each transition is %uo: the total rate of absorption of energy is 
given by multiplying these factors by the transition rate R, giving 



^ = ivr/mW(|0 nm | 2 > . (42) 

Both (^) and ( |42"1) are quadratic functions of <p e g, but it is not immediately clear how they 
can be related. We will now discuss why they are equivalent. 

The mean-square matrix element can be estimated from the correlation function C^t) 
of the effective potential 

(\<Pnm\ 2 ) = ^-Re / dt e^C^(t) (43) 
Tinn Jo 

where the correlation function is defined by 

Cab (t) = (A(r, p)B(r t , p t )) = J dv J dp A(r, p)B(r t , p t )S[E F - H(r, p)] 

Q,'(E) = Jdr J dp S[E - H(r, p)] (44) 

and r t , p t are the phase space coordinates evolved under the Hamiltonian dynamics for time 
t, starting from the point (r, p). It will be convenient to define a propagator P(r, r'; t) which 
gives the probability of reaching r' from r in time t: 
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P(r, r'; t) = (5(v t - r') )9(t) 



(45) 



where 8(t) is a step function which makes the propagator zero for negative time. The 
averaging will be defined in the next subsection. With this definition we have 

C„(t) = ^JdrJdr' P(r )r ';t)0(r)0(r') = (0|P(t)|0) (46) 



where the operator P(t) is defined by analogy with (23). Introducing the Fourier transform 
P{oj) of the propagator, we have 



\dt) = ^ V Re 1 dt ^M^M 



\ve 2 u 2 Re(0 eff |P(u;)|0 cff ) (47) 



In the next section it will be shown that there is a general relation between the propagator 
and the polarisability operator: 

U{u) = ev[T + hjP(uj)} (48) 

where v is the density of states per unit volume. A relation of this form has been given by 
Kirzhnitz [13]. We present a detailed derivation, based on Liouville's equation, below. If 
the potential V e g appearing in (|i~6|) is constant within the conducting particle, we may write 
n = uV, where V is the volume of the particle. Substituting this into fl39|) reproduces 
thus establishing its equivalence to ([4^). 



B. General relation between polarisability and the propagator 

We will now relate the polarisability operator Il(r, r'; uo) to the probability propagator 
P(r, r';t), which is the probability that an electron, released at r' with energy equal to 
the Fermi energy E-p, will be at position r after time t. The discussion will be classical; a 
quantum mechanical derivation proceeds along similar lines. 

Let the phase-space distribution be /(r, p;i): this will, when convenient, be written 
f(a, t) where a = (r, p). The Hamiltonian will be assumed to be of the form 

H(a, t) = H (a) + X^H^a) (49) 

where we will be interested in the case where H = p 2 /2m + V(r) and Hi = <f)(r). The 
perturbation parameter X{t) is assumed to be small, so that f(a,t) may be expanded as a 
series in X(t): we will be interested in the expansion as far as the first order term: 

f(a, t) = f (a) + f dt' X(t') g(a, t, f) + 0{X 2 ) . (50) 

J — oo 

Substituting into the Liouville equation, d t f = {f,H}, it is found that f is a function of 
the unperturbed Hamiltonian H (a), and that the kernel g(a,t,t') of the first order term 
satisfies 
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X(t)[g(a,t,t')-{H l J } a }+ f dt' X(t')[d t g-{H ,g}} , = (51) 

J — oo 

which is valid for all X(t). The first term implies that 

g(a,t,t) = g(a) = {Fi,./b}a = {H x , H } a ^(H (a)) . (52) 

The second term implies that dg/dt = 0, where d/dt is the total time derivative along a 
trajectory, so that 

g(ai,t,t') =g(a, t-t') = ^(H (a)) ( ^(a t „ t ,(a)) . (53) 
The required approximation is then 

f(a,t) = f (H (a)) + §§(#o(a)) f Jt' X(0^(ow) • (54) 
We will use an alternative form, obtained by integration by parts 
f{a,t) = h{H Q {a))+X{t)^{H Q {a))H x {oL) 



dfo 
dE 



(H (a)) f dt' X(t') H x {a t -t>{a)) ■ (55) 



We will assume that the integral converges. For ergodic systems this requires that the 
microcanonical average of Hi(a) vanishes. The density of available states in phase space 
is (27rh)~ d , where d is the number of degrees of freedom. For a system of fermions, the 
appropriate density function is f (a,X) = 9[E F — H(a,X)]/(27ih) d , where 9(x) is the Fermi- 
Dirac distribution, which can be approximated by a downward step function when the 
temperature is small compared to the Fermi temperature. 
Now the charge density of electrons is 

p(r, t) = e J dp f(r, p; t) . (56) 

The number density of electrons N(Ep, r) and the density of states per unit volume at the 
Fermi surface z/(i?p,r) are 

N(E F , r) = -^- d J dp 9(H (r, p) - E F ) , u(E F , r) = ^N(E Fl r) (57) 

respectively. Also, the local average of any quantity A(r, p) for electrons at the Fermi surface 
is defined as 

= ldpA( I ,p)S(H ( I ,p)-E F )/ J dp <5(ifo(r,p) - E?) 

= (2rt)A^,r) / dP ^ PWg ° (r ' P) - gF) (58) 
^From (|55"D and the definition (|56|), we have: 
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p(r, t) ~ eN(E F , r) + eX(t)u(E F , r)0(r) 

-e^ F ,r) /* dt' X(t') fdv' [dp5[E-H ( r ,p)]5[r' - rt _ t ,(r,p)](f)(r') + 0(X 2 ) 
= eN(E F , r) + eX(t)z/(E F , r)0(r) 

- eu(E F , r) T dt' X(t') [ dr' (6[r t _ t ,(r, p) - r'])^ 0(r t ,) + 0(X 2 ) (59) 

With the definition of the propagator 

P(r,r';t) = 0(t)(5[r'-r t (r,p)]} EFir (60) 
and recalling the definition of the polarisation operator, (|25|) , we find: 

n(r, r, £ - 1') = e 6(t - t')u(E F , r) 8(r - r')8(t - t') + <9 f P(r, r'; t - t') (61) 

or alternatively, in the frequency domain 

Il(r, r', u) = e is(E F , r) [8(r - r') + iuoP(r, r'; w)] . (62) 

We will introduce an operator z>, which is diagonal in the position representation, so that 
(r| v\(f)) = v(E F ,r)<f)(r). Equation fl62| ) may then be written in the form 



U{uj) =eu[T+iuP(uj)} . (63) 

There is also a relationship between the non-local conductivity £(u;) and the propagator 
P{uj), which has previously been obtained by Serota and co-workers [10,11] (with an alter- 
native derivation given in [8]). Their derivation was specific to the case of diffusive electron 
motion, whereas that given above also includes the ballistic case. 



IV. PARTICULAR FORMS FOR THE POLARIS ABILITY 

A. Spatially homogeneous, ballistic system 

For a spatially homogeneous system, II(r, r'; u) is a function of r — r', and is conveniently 
represented by its Fourier transform, II(q, uo): in d dimensions 

fi(w) = ^yj d ^ ta) n (l> ^ foil > ( r ta) = ( 64 ) 

where V is the volume of the system. In the case where the electron motion is ballistic, the 
propagator is, for d = 3, 

P(r,t';t) = ^5(R-v F t) , R = |r - r'| , 

p ^-w^ RI " (65) 

and the Fourier transform representation of the polarisability is 
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U(q,u) =ve(l- ^ lo S 



A + l 



A - 1 



which is the semiclassical limit of Eqs. (12.48a,b) in [9]. In two dimensions, U(q,uj) is given 
by 

1 + i(A - I)" 1 / 2 for A > 1 
U(q,u)=eu{ (67) 
l-fl-A 2 )- 1 / 2 for A < 1. 



B. Low and high frequency limits 

In the low frequency limit, it is immediately clear from (Q) that the induced charge 
density is p(r) = z/e[0(r) — (0)] where (0) is the space average of <fi over the particle. We 
shall be only be concerned with cases where (0) vanishes, so that we may write 

II(r, r'; uj) ~ u e 5(r - r') , u < uj c . (68) 

For sufficiently high frequencies, and sufficiently far from the boundary of the particle, the 
conductivity is local, with value a(u): 

E i:/ (r, r'; uj) = 5,^ S(r - r') o(u) (69) 

and the bulk conductivity cr(co>) may, in the case of diffusive electron motion, be approximated 
by the Drude formula 

. , ve 2 D m ON , , 

= T+l=f T= NBE D (70) 

where D is the diffusion constant, and m the electron effective mass. In the case of ballistic 
electron motion, the bulk conductivity is determined purely by the inertia of the electrons, 
and is non-dissipative: 

*{<*) = — ■ (71) 
row 



When the non-local conductivity can be approximated by (p9|) , the non-local polarisability 
takes the simple form 

n(r,r» = ^V i V:5(r-r') 

UJ 

^Mr r'). (72) 



UJ 



This approximation is expected to be valid when uj ^> u c , and when both r and r' are 
much greater than a distance A from the boundary: in the ballistic case A = vf/uj, and 



in the diffusive case A = JD/uj. The same conclusion can also be reached by considering 
the expressions flBTp, (pH) in the limit A — » 0: for d = 3 we find that H(q,u) ~ ^ue\ 



2 



vev |g 2 /3a; 2 , which is equivalent to the Fourier transform of (|72|) when the conductivity is 
given by (|7T|). 
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C. Polarisability close to a boundary 



Next we consider the polarisation charge close to the boundary of the particle. Here we 
are concerned with the high frequency case, to ^> u c . In the low frequency case ( |68D gives 
an adequate approximation, but our discussion of the high frequency case assumed that the 
conductivity could be approximated as that of a homogeneous system. Another reason for 
discussing the boundary separately is that we expect that the polarisation charge density 
may have a singularity there. 

We may assume that for uj 3> uo c the polarisability operator is short ranged. A smooth 
boundary may therefore be approximated locally by a flat surface, z = in some local Carte- 
sian coordinates. The polarisability is given by fl53|), and we approximate the propagator 



from r' = (x' , y', z') to r = (x, y, z) by the sum of a direct contribution and a contribution 
originating from an image source at r[ m = (x',y', —z'), so that 



p(r) = ve 



(r) dr> [P(r, r'; u) + P(r, r^; 



(73) 



The charge charge density is concentrated in a narrow layer at the surface, and may typically 
be approximated by writing 



p(v) = p s (z)K(S) 



(74) 



where S labels points on the boundary, and z is a coordinate normal to the boundary. In this 
case, the potential in the neighbourhood of the boundary is of the form 0(r) = <p s (z)K(S), 
where d>Jz) satisfies 



p s (z) = ve 



4> s {z) +if dz 1 [G((z - z')/K)) + G((z + z')/h)]<t> s {z') 
A Jo 



(75) 



where A = vf/uj, and the function G(x) is easily related to the Fourier transform of H(q, u). 



D. Diffusive electron motion 

In the diffusive case, it is possible to write a useful eigenfunction expansion for the 
linear response functions: for t > the propagator P(r, r'; t) satisfies the diffusion equation 
d t P = DV 2 P, or [iu — W 2 ]P(r, r'; uj) = —5(r — r'), and satisfies the Neumann boundary 
condition. It can be expressed in terms of a set of eigenfunctions Xn( r ) of the Helmholtz 
equation (V 2 + kfyxnij) = 0, satisfying the same boundary condition: n.Vxn = 0, where n 
is a normal vector on the boundary of the particle. The propagator can then be written 

^M = E ; ,.,_ 1 Dfc2 IXn)(Xn| • (76) 



n 



Expansions for other linear response functions are easily obtained in the same form. For 
example, if the density of states per unit volume v is independent of r, (51) implies that the 



polarisability can be written in this form, with the coefficient of the operator |Xn)(Xn| given 
hyveDkl/(Dkl-iu). 
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V. THE SELF-CONSISTENT FIELD 



A. Approximate equations for the self-consistent fields 

Here we discuss how the solution of the self-consistent equations can be greatly simplified 
by the use of 'semiclassical' approximations. We consider the electric dipole absorption first. 

Calculation of the electric dipole absorption coefficient via either ( j39f ) or (f47[) requires 
the self-consistent fields e fr(r, u), which is given by equation (P2|): 

|0 ext ) = [i-UTL(u)]\M . (77) 

Formally, solution of this equation requires calculation of the inverse of I — UII(uj): this 
could be done explicitly in a numerical calculation by expanding in a suitable basis set. We 
will aim instead for an approximate analytic solution. For frequencies small compared to 
the plasma frequency u p , the external electric field is 'screened' by polarisation charges, so 
that the internal field is much smaller than the externally applied field. The key physical 
intuition is that the external electric field is almost exactly cancelled by the electric field due 
to the induced charge density p(r). Let p c i(r) be the charge density induced on the particle 
by a static external field, according to classical electrodynamics: this charge density gives an 
induced electric field which precisely cancels the externally applied field inside the particle. 
For frequencies small compared to the plasma frequency, the induced charge density is well 
approximated by p c i(r): we will assume that 

\p) = IPci) + 0{uj/uj p ) + O(a/a ) (78) 

where a is the characteristic dimension of the particle, and ao is the Bohr radius. The 



classical charge distribution formally satisfies an equation analogous to (|77|), in which the 
term representing the internal field \4> e s) is set equal to zero: 

|0ext) + U\Pcl) = • (79) 

We will denote our approximation to the effective potential |0 e fr) by |0): it is the potential 
which generates the polarisation charge \p c \), and is given by 

\p cl ) = fl(u)\<f ) ) (80) 

or equivalently by |0 cxt ) + Ufl(u)\(f)) = 0. Comparing with ([77Q, it is clear that this solution 
\(f>) is a good approximation to \4> e s) provided ||LTI(u;)|| 1, where ||X|| is an appropriate 
norm of the operator X. To estimate this norm, we consider the effect of an arbitrary 
potential 0: at zero frequency, the induced charge density is 110 = ez/0, and for a particle 
of characteristic dimension a in d dimensions, the induced charge may be approximated 
by a dipole formed by charges of magnitude Q ~ pa d , with separation a: this results in 
an electrical potential of magnitude 0' ~ eQ/(eoa). In three dimensions this leads to the 
following estimate for ||f7II|| ~ 0'/0 

UU (81) 
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where u v is the three dimensional bulk plasma frequency, u p = [Ne 2 / (AireQin)] 1 ^ 2 , N being 
the electron density. For frequencies uj 3> tv c , a similar argument gives 

uj 2 

\\UU(u »w c )|| . (82) 

UlUlr 



In the case of ballistic dynamics uj c = vp/a, equations (81) and (|82D are of the order of 
cr /(e a; c ) and a /(e uj) respectively. 

We can, in principle, determine improved approximations to the exact solution of (|78|) 
from the solution of (p0|): 

\M = -[i-ufi(u)]- 1 ufi(u)\<t>) 

= \<t>)-[i-UU{u)]- x \<i>) . (83) 

Equations ( |5TD and fl82|) show that the correction term in (|83|) is small. Having shown that 
|0eff) ~ |0) ; we discuss how to estimate the solution |0) of ( jSOf) in Sees. VB and VC. 

Finally we consider the semiclassical solution for the magnetically induced field, E ind 
which satisfies flBT]): 

V.[/-ia;eoE(a;)]|E ind ) = 0. (84) 

At low frequencies we can estimate the conductivity by E ~ Ne 2 /mu s where u s is the 
scattering frequency; at high frequencies u is replaced by the frequency uj. In the low 
frequency limit we therefore estimate 

— i|EH||~^. (85) 

In the frequency range that we are concerned with, the term involving the identity operator 
in (|84]) is therefore negligible; the same conclusion applies when uj ^> u s . We can therefore 
find an approximate solution to fl8"4| ) by requiring that the induced charge density is zero, 
i.e. 

V.EM|E ind ) = . (86) 

This justifies the neglect of a cross-term in (||) with an electric dipole induced by the magnetic 
field. A solution to (^) may be determined by choosing an initial approximation E' nd which 
satisfies VAE[ nd = io;B ext . A polarisation charge p[ nd would be generated from this field. An 
additional field which is the gradient of a potential x( r ) is added, such that E ind = E[ nd + V%- 
The condition upon x f° r (0) to be satisfied is 



\pU = fi(u;) |x) • (87) 



This equation is analogous to 
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B. Solutions in the ballistic case 



We will discuss approximate solutions of (^0) valid in the limits w<Cw c and u ^> u c . The 
first of these represents the static potential required to hold the classical charge distribution 
in place in the zero- frequency limit: it will be written |0 s t a t)> an d its form is immediately 
apparent from (|63|): 

l^stat) = — |Pcl) • (88) 

This is simply a linearised Thomas- Fermi approximation [12]. The semiclassical approxima- 
tions underlying this expression assume that the potentials are slowly varying on the scale 
of the Fermi wavelength. This is questionable when the charge density has a singularity on 
the surface of the particle: we will return to this later. 

In the limit u) 3> oo c , we found [equation (0)] that the polarisability may be approximated 
by II ~ — (\a{uj) / uj)V18{t — r'), for r and r' not too close to the boundary. For points not 
too close to the boundary, or where the charge density is non-analytic, we can approximate 
the solution of fl80|) by a 'dynamic' potential, which is of the form |0d yn ) = A|-?/>), where 
V 2 ^(r) = p c i(r), and n.V^r) = 0. Substituting these forms into (p0|) we find that A = 
—\uj/a{uj): 

|0dyn)=-^TiV) , VV) = | Pcl ). (89) 

An interpretation of the dynamic potential is that it moves the polarisation charge into 
place. 

Close to the boundary (159) is not necessarily a good approximation to the effective 



potential. One reason is that the approximations underlying ([72]) fail, and the polarisation 
must be described by Q73|) or (]75|). Another reason is that the charge density \p c \) has a 
singularity there. In the notation of (|74"D, the projected charge density p s (z) is concentrated 
on the boundary in the three dimensional case [4] so that p s {z) ~ ^{ z )i and in the two 
dimensional case it diverges on the boundary, such that p s (z) ~ z~ x l 2 for z > [6]. The 
form of (|75|) indicates that the potential (f) s {z) also has the same type of singularity as the 
charge density at the boundary. 

To summarise, the following picture emerges. For low frequencies, u <C uj c the potential 
is approximately sta t(r). At high frequencies, u ^> u> c , the potential is well approximated 
by 0d y n(r) within the interior of the particle. In the vicinity of the boundary, the potential 
has a dominant divergent contribution, which is well approximated by s tat(r)- 

C. Diffusive case 



In the diffusive case, equation (p0|) can be solved exactly, using the representation of the 
propagator in the form (|7B|). Expanding the potential \<p) in terms of the functions \x n ) 
leads to the expression 



1 



J2(Xn\Pcl) 1 - 7^ \Xn) = |&tat) + I0dyn) (90) 

< 

dynamic contributions [7]. 



Note that in the diffusive case the potential is precisely equal to the sum of the static and 
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VI. DISCUSSION: CALCULATION OF THE ABSORPTION COEFFICIENT 



Once an adequate approximation for the effective potential has been obtained, the elec- 
tric absorption coefficient is obtained from P7|) : the absorption coefficient is proportional 
to u 2 (<j)\P(uj)\(j)). Previous papers [4-8] have discussed methods for the evaluation of the 
absorption coefficient using equation (|47|) in some detail, for specific cases. This paper has 
presented a general approach to the determining of the effective potential, and some remarks 
on applying this to calculating the absorption coefficient may be useful. 

When the electron motion is diffusive, the absorption is very easily evaluated using fl90|) 
and QTSp. It is found that the coefficient is proportional to u 2 , and that (at least within the 
framework of the approximations used in section V) the frequency scale uj c plays no role. 
The absorption coefficient can be shown to be exactly equal to the classical value in this 
case [7,5]. 

The case of ballistic electron motion is more difficult. It might be expected that C^iui) = 
{(j)\P{uj)\4>) approaches a non-zero limit as uj — > 0, implying that the absorption coefficient 
is proportional to uj 2 for low frequencies. This expectation is correct for cases where the 
electron motion is ergodic (the most important cases being diffusive electron motion, and 
the ballistic case with a rough boundary). In the case of integrable electron motion, which 
can be realised experimentally if the boundary appears smooth on the scale of the Fermi 
wavelength, C^u) typically approaches zero as uj — > in a non-analytic fashion. In the 
important special case of particles with circular symmetry, C^iuj) is zero for uj < u> c , 
where uj c is the frequency of a glancing circular orbit. Thus the electron dynamics plays an 
important role in determining the low frequency absorption. 

The high frequency absorption, by contrast, is determined by the nature of the singu- 
larities of the function f(t) = 0(r t ), which can result from singularities in the motion r(t), 
or in the potential 0(r). We will discuss the ballistic case. The dominant contribution 
comes from the singularities of </>(r) in the neighbourhood of the boundary. The absorption 
coefficient was calculated in [4] and [6] for the three dimensional and two dimensional cases 
respectively, assuming that the effective potential is equal to sta t( r )- In three dimensions, 
this potential is a delta function singularity concentrated on the boundary, and in two di- 
mensions it diverges as z~ 1//2 , where z is the normal distance from the boundary. These 
forms for the potential imply that (<p\P(uj)\<j)) ~ uj° in three dimensions and ~ uj' 1 in two 
dimensions. In [4] and [6] it was predicted that in the case of particles with circular sym- 
metry, the absorption coefficient shows a sequence of resonances superposed on a regular 
contribution increasing as uj 2 and uj for three and two dimensions respectively. 

The more sophisticated approach introduced in this paper indicates that the potential 
|0) should satisfy (|30D, whereas the potential used in [4], [6] was simply the static potential, 
satisfying \p c \) = n(O)|0 st at)- We must discuss the extent to which this refinement will 
change the results. In section V B we argued that the singularities of the effective potential 
at the boundary are the same as those of the 'static' potential. We can therefore hypothesise 
that the more refined theory would make a quantitative rather than qualitative difference 
to the results. This issue will be addressed in a subsequent paper, which will consider the 
inversion of (|75|) to determine the potential <p s (z) from the charge density p s {z), and its use 
to estimate the high frequency absorption coefficient. 
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